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Abstract. The spectral libraries of hot, massive stars which are 
implemented in the population synthesis code Starburst99 are dis- 
cussed. Hot stars pose particular challenges for generating libraries. 
*~* , They are rare, they have an intense radiation field and strong stellar 

winds, and a luminosity bias towards ultraviolet wavelengths. These 
properties require the utilization of theoretical libraries. Starburst99 
uses static non-LTE models at 0.3 A resolution in the optical, spheric- 
ally extended, expanding models at 0.4 A resolution, in the satellite- 
ultraviolet, and blanketed, low-resolution radiation-hydrodynamical 
models in the extreme ultraviolet down to X-rays. I review the main 
features of each library, compare them to observations, and discuss 
their link with stellar evolution models. 



Keywords : stars: atmospheres - stars: early-type - stars: winds, 
outflows - galaxies: stellar content 



Introduction 



The past decade has witnessed considerable progress in charting the demo- 
; "*! ] graphics of high-redshift galaxies. Panchromatic surveys have established the 

luminosity functions of star-forming galaxies both in the ultraviolet (UV) 
(Reddy & Stcidcl 2009) and in the infrared (IR) and sub-mm (Wardlow et al. 
2011). Intermediate dispersion spectroscopy of bright Lyman-break galaxies 
(LBGs) has been pivotal in stellar population trends that cannot be identified 
from photometric data alone. Such data can now be routinely obtained for 
"normal" galaxies at redshifts as high as z f=s 4, which corresponds to an age 
of the universe of ^900 Myr (Shapley 2011). Owing to the young ages, the 
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rcstframc-UV spectra of these galaxies — redshifted into the optical — show 
the signatures of hot, massive stars. The restframe-UV spectra of LBGs exhib- 
its striking resemblances to those of star-forming galaxies in the local universe 
(Schwartz et al. 2006). Both distant and nearby star- forming galaxies and 
their hot-star populations are the targets of the spectral libraries discussed 
in this paper. The stellar contribution to the spectral energy distribution of 
these galaxies falls into three categories: (i) the restframe-UV with its distinct 
stellar-wind absorption lines from O stars; (ii) the optical/near-IR with nu- 
merous weak photosphcric absorptions from B and A stars, plus conspicuous 
emission lines from H II regions; (iii) the ionizing extreme UV, which is dir- 
ectly related to the UV wind lines and the optical nebular lines. I will use this 
astrophysical categorization of the spectra of star-forming galaxies to discuss 
the groups of spectral libraries of hot stars as implemented in the evolutionary 
synthesis code Starburst99 (Leitherer ct al. 1999; Vazquez & Leitherer 2005; 
Leitherer & Chen 2009). 



2. Why use Theoretical Libraries? 

Spectral libraries of hot, massive stars must be built from theoretical model 
atmospheres. This fact is driven not by choice but by nature. Massive stars 
are rare and therefore distant — the closest early-0 star £ Puppis is at a dis- 
tance of ^450 pc. This, together with the concentration of massive stars in the 
Galactic plane, leads to large dust extinction in the optical/UV and associated 
large flux uncertainties. Moreover, large gas column densities produce strong 
interstellar absorption lines which can blend stellar features or even fully de- 
press the stellar continuum in the far-UV. Pellcrin ct al. (2002) used the Far 
Ultraviolet Spectroscopic Explorer (FUSE) satellite to collect far-UV spectra 
of OB stars in the Galaxy, the LMC, and the SMC. Galactic stars even with 
moderate dust reddening as low as E{B — V) = 0.1 have H2 columns which 
make large portions of the spectra unusable for library purposes. Of course, the 
foremost reason for having to rely on theoretical libraries is the lack of metal- 
poor, hot, massive stars with a metal abundance of Z < jkZq in the local 
universe. Massive stars are short-lived and are forming during an epoch when 
local star-forming galaxies (including the Milky Way) have evolved from low to 
high metallicity. All known metal-poor star-forming galaxies are at distances 
beyond ^10 Mpc, which rules out observations of individual massive stars for 
the purpose of generating libraries: not only are the stars dauntingly faint, 
but they are located in dense star clusters whose members cannot be resolved 
even by spectrographs in space. This dictates a hybrid approach when build- 
ing libraries of hot massive stars. We must largely rely on model atmospheres 
but whenever possible, cross-checks and comparisons with observations must 
be performed to validate this approach. 

Modeling the atmospheres of hot stars is challenging due to severe de- 
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partures from local thermodynamic equilibrium (LTE) caused by the intense 
radiation, low densities, and the presence of supersonic stellar winds initiated 
by the transfer of momentum from the stellar radiation field to the atmo- 
spheric plasma. Fortunately, rapid progress during recent years has led to an 
astounding degree of sophistication in the latest generation of models (e.g., Puis 
2008). The main challenges are in the areas of non-LTE, line-blanketing, and 
radiatively driven winds. Several groups have independently developed model 
atmospheres for OB stars, with a different emphasis on one or more of these as- 
pects. Plane-parallel models, such as ATLAS (Kurucz 2005) or TLUSTY (Lanz 
& Hubeny 2003, 2007), cannot account for spectral lines forming in the wind 
but are suitable for photospheric modeling. The major spherical wind model 
atmospheres are PoW-R (Hamann & Grafcner 2004), PHOENIX (Hauschildt 
et al. 1999), CMFGEN (Hillier & Miller 1998), WM-Basic (Pauldrach et al. 
2001), and FASTWIND (Puis et al. 2005). PoW-R and PHOENIX are optim- 
ized for the modeling of Wolf-Rayet (W-R) and late-type stars, respectively. 
FASTWIND 's main application is the computation of optical and near- infrared 
H and He lines. On the other hand, CMFGEN and WM-Basic are both widely 
used for modeling the UV spectra of hot stars. 

Errors due to missing or incorrect atomic data are a concern for the compu- 
tation of spectra (Hillier 2011). However, this is less of an issue in hot stars as 
compared to their cool counterparts where the presence of dust and molecules 
adds to the complexity of the modeling. Overall, the state of spectral model- 
ing of hot stars gives reasonable confidence in the reliability of the computed 
spectral libraries. 



3. The Ultraviolet Spectral Region 

The space-UV part of the spectrum (defined as the wavelength region between 
912 and 3000 A) is dominated by strong stellar-wind lines from hot massive 
stars. Massive stars are luminous and develop strong mass loss already when 
on the main-sequence. The resulting stellar winds are radiatively driven, with 
radiative momentum being transferred into kinetic momentum via absorption 
in metal lines (Puis, Vink, & Najarro 2008). Since there is a well-defined 
relation between the radiative and the kinetic momentum (the efficiency rate 
is about 30%), the kinematic properties of the wind (as observed in the wind 
lines) reflect the basic stellar properties, such as effective temperature T e s, 
luminosity L, and gravity \ogg. 

We used WM-Basic for modeling the UV spectra of hot stars (Leitherer 
et al. 2010). A complete model atmosphere calculation consists of three main 
blocks: (i) the solution of the hydrodynamics; (ii) the solution of the non-LTE- 
modcl (calculation of the radiation field and the occupation numbers); (iii) the 
computation of the synthetic spectrum. The three cycles are interdependent 
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Figure 1. Examples of theoretical UV spectra. Left: T c g = 37,800 K; logg = 3.7; 
Z = Zq. The corresponding spectral type is 06 III. Right: T c g = 15,300 K; logg = 
1.9; Z = Zq. The corresponding spectral type is B3 la. The most important spectral 
lines are identified in both spectra. From Leitherer et al. (2010). 



and must therefore be solved iteratively. In the first step, the hydrodynamics 
is solved for a set of T c g, logg, radius R, and Z, together with an assumed 
radiative line acceleration. In a second step the hydrodynamics is solved by 
iterating the complete line acceleration (which includes the opacities of all 
important ions) and the temperature and density structure. Afterwards, these 
structures are again iterated together with the line acceleration obtained from 
the spherical non-LTE model. The final step consists of the computation of 
the synthetic spectrum. 

Two examples are shown in Fig. [TJ The hot star (left panel) has a spectral 
type of 06 III. Prominent spectral lines are identified. Most of these lines 
are formed in the stellar wind, as can be seen from the blueshifted absorption 
components with velocities exceeding ^2000 km s _1 . S VI A939, O VI A1035, 
P V A1123, N V A1240, and C IV A1550 are the strongest features that are uni- 
formly present in O stars. Si IV A1400 is weak in this particular model but can 
become a strong line in supergiants whose denser winds lead to recombination 
from Si 4+ (which is the dominant ionization stage) to Si 3+ (Walborn & Panek 
1984; Drew 1989). The wavelength region longward of 1800 A is devoid of both 
stellar- wind and photospheric lines. This region has little diagnostic value for 
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constraining an O-star population. The O-star spectrum in Fig. 1 can be con- 
trasted with the B-star spectrum in the right panel which corresponds to a 
B3 la supergiant. The high-ionization lines present in the O-star spectrum are 
much weaker or completely absent in B stars. The strongest lines are, among 
others, C III A1176, C II A1335, Si IV A1400, and Fe III A1893. This star has 
a terminal wind velocity of v^ = 300 km s _1 , resulting in comparatively small 
blueshifts of the absorption components. The line-blanketing by photospheric 
absorption lines is significant, in particular at the shortest wavelengths where 
the rectified flux is close to the zero level. 

Using this method, we generated a spectral library of 430 spectra for in- 
clusion in evolutionary synthesis models of star-forming galaxies. The chosen 
stellar parameters cover the upper left Hertzsprung-Russell diagram at L > 
1O 275 L and T e g > 20, 000 K. The adopted elemental abundances are 0.05 Zq, 
0.2 Zq, 0.4 Zq, Zq, and 2 Zq. The spectra cover the wavelength range from 
900 to 3000A and have a resolution of 0.4 A. We compared the theoretical spec- 
tra to data of individual hot stars in the Galaxy and the Magellanic Clouds 
obtained with the International Ultraviolet Explorer (IUE) and FUSE satel- 
lites and found very good agreement. We then implemented the library into 
the synthesis code Starburst99 where it complements and extends the empirical 
libraries. 

We generated a model series for a standard single stellar population, i.e., an 
instantaneous burst of star formation following a Salpeter initial mass function 
(IMF) with mass limits of 1 and 100 M Q and following standard evolutionary 
tracks. In Fig. [2j we show the comparison with the IUE library at Zq. The 
spectral resolution of the WM-Basic library was decreased by a factor of 2 to 
match that of the IUE library. Plotted are spectra covering ages from 1 to 
20 Myr at 9 time steps. 20 Myr is the approximate evolutionary time scale of 
a 10 Mq star. Both spectral libraries are complete at (and beyond) this age. 
Evidently the overall agreement between the spectra generated with the two 
sets of libraries is excellent for all ages shown in this figure. The age evolution 
of the population is most clearly reflected in the weakening of the N V A1240 
and C IV A 1550 P Cygni profiles with time, which results from the decreasing 
O-star contribution. Si IV A 1400 exhibits the well-known luminosity effect and 
is strongest when massive O supergiants appear between 3 and 5 Myr. N IV 
A1720 behaves similarly. Numerous narrow absorption lines in the empirical 
spectra have interstellar origin. Obviously such contamination is absent in 
the theoretical spectra, thereby greatly facilitating automatic spectral fitting 
procedures that utilize the full spectrum. Notice Si IV, whose interstellar 
components contribute significantly to the total profile and introduce a major 
uncertainty in the interpretation of this line. Lyman-a has strong damping 
wings in the empirical data due to Galactic H I absorption. This affects the 
blue edge of the N V P Cygni profile. Again, the library generated with 
WM-Basic does not have this issue. Finally, the normalization of the empirical 
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Figure 2. Comparison of synthetic spectra of a single stellar population obtained 
with the WM-Basic library with Z = Zq (thick lines) and an empirical IUE library 
of Galactic stars (thin) between 1200 and 1800 A. Each spectrum represents one age 
step from 1 to 20 Myr as labeled on the right. Salpeter IMF with mass limits 1 and 
100 Mq. From Leitherer et al. (2010). 

spectra around 1450 A is questionable. Severe blanketing makes the continuum 
determination in this wavelength region a challenge and causes the incorrect 
continuum location in the IUE data. 

The new theoretical library is a significant improvement over the existing 
empirical library and allows the modeling and interpretation of the UV spectra 
of massive star populations for essentially any chemical composition of interest 
and at any redshift. 



4. The Optical Spectral Region 



The WM-Basic library as implemented into Starburst99 is confined to wavelengths 
below 3000 A. The WM-Basic atmosphere code makes simplifying assumptions 
which make it not suitable for calculating the optical line spectra (such as the 
omission of Stark broadening, which is important for computing the Balmer 
lines) . The optical spectral region of a star- forming galaxy is weighted towards 
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Figure 3. Comparison between non-LTE and LTE models for solar abundances. 
Dashed: Kurucz (LTE) at a low resolution of 20 A; dotted: Kurucz (LTE) recomputed 
with SYNSPEC at a resolution of 0.3 A; solid: TLUSTY (non-LTE) and SYNSPEC 
at a resolution of 0.3 A. The non-LTE models are a dramatic improvement over 
previous model spectra for hotter T e g. From Martins et al. (2005b). 



less massive stars, which allows alternative models to be used. Whereas the 
UV wind lines are related to O stars with masses of 20 to 100 M@, the op- 
tical continuum and line spectrum traces less massive stars of spectral types 
B and A. Wind effects are less pronounced in these stars, but non-LTE and 
line-blanketing are crucial. Consistently treating metal line-blanketing in non- 
LTE atmospheres is challenging but feasible. Line-blanketing has the effect of 
blocking UV radiation, which then emerges at optical and longer wavelengths. 
Some of the blocked photons are back-scattered, causing a back-warming ef- 
fect, which leads to an enhanced ionization in the inner atmosphere. This 
important process has now been incorporated into the most recent generation 
of models and has led to the realization of a significantly revised T e g scale of 
OB stars (Martins, Schaerer, & Hillier 2005a). 

Martins et al. (2005b) computed 1654 high-resolution stellar spectra with 
a sampling of 0.3 A and covering the wavelength range from 3000 to 7000 A. 
The grid covers the full Hertzsprung-Russell diagram for chemical abundances 
of 0.1 Zq, 0.5 Zq, Zq, and 2 Zq. The spectra were generated using the most 
appropriate models available for each temperature and gravity value. The lib- 
rary uses non-LTE line-blanketed models for hot stars (T e g > 27, 500 K) and 
PHOENIX LTE line-blanketed models for cool stars (3000 K < T cS < 4500 K). 
Kurucz atmosphere models were used for the parameter space in between. The 
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replacement of the traditional Kurucz atmospheres at the hot and cool end 
is significant. PHOENIX models account for the overall energy distribution 
in a self-consistent way because of their extensive molecular line list. These 
molecular transitions are missing in the widely used Kurucz models. Lejeune, 
Cuisinier, & Buser (1998) provided an empirical correction to the Kurucz mod- 
els. This correction becomes obsolete with the new fully blanketed models. 
Important age diagnostics such as the helium lines are strongly affected by 
non-LTE effects in OB stars. LTE models predict too weak equivalent widths. 
An example of the corresponding improvement is shown in Fig. [3] The fig- 
ure suggests that for higher temperatures non-LTE models make a significant 
difference, as already shown by Hauschildt et al. (1999). LTE models are 
adequate for B-type stars, whereas for hotter stars non-LTE effects become 
progressively more important. 

Gonzalez Delgado et al. (2005) tested this library by linking them to 
a population synthesis code to calculate high-resolution spectra of a stellar 
population in the 3000 - 7000 A range. These models predict the indexes of 
metallic lines and the absorption-line profiles of the hydrogen Balmer series and 
helium lines as a function of age and metallicity for an instantaneous burst. 
The resulting synthetic models agree very well with observations of individual 
star clusters. Discrepancies are mainly associated with the limitations and 
uncertainties in stellar evolution models. These limitations can be relevant 
when red supergiants become important in their contribution, in particular at 
low metallicity. 



5. The Ionizing Ultraviolet Region 

The wavelength range below 912 A is the source of the radiative heating of the 
interstellar gas. While this wavelength region is usually not considered part 
of a spectral library, it is relevant in that it provides an additional invaluable 
diagnostic. The stellar radiation field in the extreme UV powers the stellar 
winds and is responsible for the strengths and radial velocities of the UV wind 
lines. Therefore the nebular emission lines in the optical spectrum (which we 
use as a proxy for the ionizing radiation field) probe the same energies as the 
wind lines and should lead to consistent results. 

The dramatic effects of non-LTE and sphericity effects in the model atmo- 
spheres are illustrated in Fig. @] where WM-Basic models are contrasted with 
classical, static Kurucz LTE atmospheres and with CoStar models (Schaerer & 
de Koter 1997), which are similar to WM-Basic models, except for the exclu- 
sion of line-blanketing. WM-basic is considered to have the most appropriate 
physical ingredients among the three models shown and is a viable choice for 
photoionization modeling. All models agree rather well in the wavelength range 
where the hydrogen ionizing radiation is emitted (A > 504 A). The number of 
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Figure 4. Comparison of the ionizing UV fluxes for an early-0 supergiant from a 
WM-Basic (solid), a CoStar (dashed), and a Kurucz model (dotted) at Z — 0.2Zq. 
From Smith et al. (2002). 



photons predicted to be emitted in the hydrogen Lyman continuum is a robust 
quantity which has changed by less than 0.3 dex since the classical work of 
Panagia (1973). The behavior of the ionized helium continuum below 228 A 
is in sharp contrast to that of hydrogen. Now, wind effects drastically alter 
the spectrum, depending on blanketing and wind density. For the parameters 
chosen in Fig. |H a vast flux excess over Kurucz models is predicted for WM- 
Basic models, yet blanketing somewhat lowers the output below that of the 
CoStar models. The behavior of the neutral helium continuum at wavelengths 
between 228 and 504 A is bracketed by the cases of H and He + . The astrophys- 
ical reason for the flux increase of the extended models is the velocity gradient 
due to the outflowing matter. The resulting Doppler shift favors less opacity 
and more free-bound emission. 

Dopita et al. (2006) and Levesque, Kewley, & Larson (2010) used the Map- 
pings III photoionization models to produce theoretical emission-line spectra of 
H II regions and star-forming galaxies with WM-Basic SEDs as input sources. 
Their extensive suite of models is intended to facilitate detailed studies of star- 
forming galaxies and their interstellar medium properties, but also serves as an 
important test of the model atmospheres. The model sets of the two groups of 
authors are equivalent but differ in some important aspects. The Dopita et al. 
models take the ratio of the mass of the central aging star cluster to the gas 
pressure as the free parameter. Therefore the ionization parameters becomes 
a function of age, as appropriate for individual H II regions. To model the 
spectra of star- forming galaxies as done by Levesque et al., the individual H II 
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regions are integrated for each model age, essentially considering star-forming 
galaxies to have spectra that consist of contributions from multiple H II regions 
at different ages. 

Levesque et al. (2010) examined the agreement between the photoioniz- 
ation models and the spectra of local star-forming galaxy populations from 
several large data sets, including the Sloan Digital Sky Survey. In general the 
models do a reasonably good job in reproducing most metallicity-sensitive line 
ratios. Some lines, however, require a harder radiation field than predicted 
by the models. This suggests that further systematic changes are required in 
the stellar population synthesis models to produce the far-UV ionizing spectra. 
It is not clear whether the model atmospheres or some other ingredient is to 
blame. Currently it appears that the adoption of stellar evolutionary tracks 
which include the effects of rotation may provide the necessary hardening of 
the spectrum (Leitherer & Ekstrdm 2012). 



6. Applications 



As a demonstration of the library's capabilities I am giving two recent applic- 
ations, one related to the chemical composition of high-redshift galaxies, and 
the other to the IMF of local star-forming galaxies. 

A major challenge in deriving abundances from nebular emission lines is 
that the strongest and best calibrated features are all located in the rest- 
frame optical, which is redshifted into the near-IR at cosmological redshifts. 
Observations in this wavelength regime from the ground are plagued by high 
(and variable) sky background and by the lack of multi-object spectrographs. 
There is therefore a strong incentive to test abundance measures based on the 
spectral features from massive stars which dominate the rest-frame UV spectra 
of high-z galaxies. 

Initial applications of this technique turn out to be quite promising (Rix 
et al. 2004; Halliday et al. 2008). A spectrum created by adding together 
75 star-forming galaxy spectra at z « 2 from the GMASS survey was used to 
estimate the iron abundance by measuring the strengths of weak photospheric 
absorption lines. A Starburst99 model with a stellar chemical composition 
of ^20% solar was found to match the data. This abundance is similar to 
those derived in individual LBGs, such as the proto-type MS1512-cB58 (a 
strongly lensed LBG at z — 2.73), for which a stellar metallicity of 0.4 Z© was 
found. The latter value is in excellent agreement with the metallicity of the 
interstellar medium of this galaxy. This method can provide an additional way 
to determine the metallicity of high-z galaxies, particularly at redshifts z > 4 
where emission-line diagnostics are difficult to apply from the ground. 
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Figure 5. Comparison between the average spectrum of 46 local star- forming 
regions generated by Leitherer et al. (2011; thick line) and a synthetic spectrum 
with a Kroupa/Salpeter IMF populated with stars up to M — 100 Mq (thin line). 
The narrow absorption lines (e.g., at 1260 A, 1335 A, 2600 A, 2800 A) are interstellar 
and absent in the synthetic spectrum. There is strong geocoronal and Milky Way 
Lyman-a contamination in the data. Note the excellent agreement for the wind 
features of N V A1240 and C IV A1550. 



The UV lines are useful for probing the IMF. The lines depend on the 
wind properties, which in turn depend on the stellar luminosity L. Since there 
is a well established mass-luminosity relation in stars, the wind lines trace 
stellar mass and ultimately the IMF of the population (Leitherer 2011). This is 
equivalent to the nebular Ha method for determining the IMF but circumvents 
the interstellar radiative transfer in gas and dust because the wind absorption 
occurs in the outskirts of the star itself. The only free parameters in the 
modeling are the star- formation history of the population and the IMF. Since 
massive stars have evolutionary timescales of less than 0.1 Gyr, star- formation 
equilibrium is reached quickly, and the population properties become time- 
independent in statistically large systems. This essentially leaves the IMF as 
the only free parameter. 



Fig. [5] is a comparison between an average UV spectrum of local star- 
forming galaxies and a standard model using the theoretical UV library. The 
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observational high-S/N spectrum is the mean of 46 individual spectra. The 
average oxygen abundance of log(0/H) + 12 = 8.5 is close to, but slightly 
below the solar value. The synthetic UV spectrum was generated for a standard 
equilibrium population at solar chemical composition following a Kroupa IMF 
between 0.1 and 100 M . (Note that this IMF is identical to the Salpeter IMF 
for massive stars.) The key constraints are the N V, Si IV, and C IV lines. 
The interpretation of the Si IV doublet requires care because of the presence 
of additional strong interstellar lines which are obviously absent in the model. 
IMF slopes steeper than ~2.6 and flatter than ^2.0 can be excluded. There 
are no additional adjustable parameters. 



The Road Ahead 



Theoretical model atmospheres of hot stars have reached a stage of maturity 
to allow the generation of astrophysically useful stellar library spectra. While 
comparison with data gives confidence in the models, I would like to end this 
contribution by mentioning known deficiencies and hurdles. 

A major issue plaguing the computed UV lines concerns wind inhomogen- 
eities. The models discussed here assume smooth winds, yet the ouflows are 
known to have significant density fluctuations (Puis et al. 2008). This causes 
some lines, in particular those from ions far from the dominant ionization stage, 
to be discrepant. An example is the O V A1371 line, which is predicted to be 
present in very hot O stars but in reality is hardly ever observed. The reason 
is understood: on average, a larger mass fraction of the winds is confined to 
higher density than assumed in the models (which are clump- free). This shifts 
the average ionization balance to lower stages and drastically weakens O V 
(Bouret et al. 2003). Efforts to account for clumped winds are underway 
(Sundqvist et al. 2011) and will hopefully allow us to address one of the most 
important deficiencies of the current spectral library. 

Completeness is an issue for almost every spectral library. In the case 
of hot stars, the challenges are both at the low- and at the high-mass end. 
At low masses, it is important to realize that evolved hot stars are not in- 
cluded. While such stars are not relevant in the presence of hot massive stars, 
they may become a significant contributor to the blue/UV luminosity once 
massive stars disappear in an ageing stellar population (Flores-Fajardo et al. 
2011). Theoretical efforts are still lagging in accounting for these stars. At 
high masses, we have to deal with Wolf-Rayet stars whose peculiar evolution- 
ary states, powerful stellar winds, and anomalous chemical abundances pose 
formidable challenges. Wolf-Rayet stars, while luminous, are often believed to 
be too few in number to be significant for inclusion in spectral libraries. How- 
ever, Crowther et al. (2010) demonstrated that the Wolf-Rayet-like stars in 
the center of the 30 Doradus Nebula in the LMC have masses possibly reaching 
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up to 300 Mq. Therefore these stars make a disproportionate contribution to 
the global UV-light budget and must be accounted for in spectral libraries. 
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